home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 18 / CU Amiga Magazine's Super CD-ROM 18 (1997)(EMAP Images)(GB)[!][issue 1998-01].iso / CUCD / Utilities / Laplace / Projects / physic_PI_24 < prev    next >
Text File  |  1996-06-28  |  2KB  |  60 lines
  1. #LAPLACE
  2. #USEFLOAT=FALSE
  3. #SIMPLIFY=TRUE
  4. #ERRBND=20
  5. #EXPDSP=7
  6. #PREC=2
  7. #F1=),R_a(
  8. #F2=g(x)=
  9. #F3=h(x)=f(x)/g(x)
  10. #F4=matrix(
  11. #F5=vector(
  12. #F6=diff(f)
  13. #F7=diff(g)
  14. #F8=diff(h)
  15. #F9=taylor(f,0)
  16. #F10=taylor(f,1)
  17. #F11=
  18. #F12=
  19. #F13=
  20. #F14=
  21. #F15=
  22. #F16=
  23. #F17=
  24. #F18=
  25. #F19=
  26. #F20=
  27. R_VM=15000,R_AM=4000
  28. dI=0.015*50e-6,dU=0.15
  29. R_a(U,I)=U/(I-U/R_VM)
  30. R_b(U,I)=(U-R_AM*I)/I
  31. dR_a(U,I)=abs(diff(R_a,I))*dI+abs(diff(R_a,U))*dU
  32. dR_b(U,I)=abs(diff(R_b,I))*dI+abs(diff(R_b,U))*dU
  33. R_a(6,16.5e-3)
  34. dR_a(6,16.5e-3)
  35. R_b(8,1.2e-6)
  36. dR_b(8,1.2e-6)
  37. R=set(R_a(10,0.095),R_a(9,0.090),R_a(8,0.084),R_a(7,0.078),R_a(6,0.072),R_a(5,0.064),R_a(4,0.057),R_a(3,0.048),R_a(2,0.040),R_a(1,0.028),R_a(0.5,0.021),R_a(0.2,0.014))
  38. dR=set(dR_a(10,0.095),dR_a(9,0.090),dR_a(8,0.084),dR_a(7,0.078),dR_a(6,0.072),dR_a(5,0.064),dR_a(4,0.057),dR_a(3,0.048),dR_a(2,0.040),dR_a(1,0.028),dR_a(0.5,0.021),dR_a(0.2,0.014))
  39.  
  40. dR_1=set(dR_a(11.5,0.091),dR_a(11,0.087),dR_a(10,0.080),dR_a(9,0.072),dR_a(8,0.063),dR_a(7,0.056),dR_a(6,0.048),dR_a(5,0.040),dR_a(4,0.032),dR_a(3,0.024),dR_a(2,0.016),dR_a(1,0.009))
  41. R_1=set(R_a(11.5,0.091),R_a(11,0.087),R_a(10,0.080),R_a(9,0.072),R_a(8,0.063),R_a(7,0.056),R_a(6,0.048),R_a(5,0.040),R_a(4,0.032),R_a(3,0.024),R_a(2,0.016),R_a(1,0.009))
  42. average(R_1)
  43. eval(error(R_1))
  44.  
  45. dR_2=set(dR_b(11.5,0.087),dR_b(11,0.084),dR_b(10,0.077),dR_b(9,0.069),dR_b(8,0.061),dR_b(7,0.054),dR_b(6,0.046),dR_b(5,0.038),dR_b(4,0.031),dR_b(3,0.024),dR_b(2,0.016),dR_b(1,0.008))
  46. R_2=set(R_b(11.5,0.087),R_b(11,0.084),R_b(10,0.077),R_b(9,0.069),R_b(8,0.061),R_b(7,0.054),R_b(6,0.046),R_b(5,0.038),R_b(4,0.031),R_b(3,0.024),R_b(2,0.016),R_b(1,0.008))
  47. average(R_2)
  48. eval(error(R_2))
  49.  
  50. dR_2a=set(dR_a(11.5,0.82e-3),dR_a(11,0.79e-3),dR_a(10,0.72e-3),dR_a(9,0.64e-3),dR_a(8,0.58e-3),dR_a(7,0.50e-3),dR_a(6,0.44e-3),dR_a(5,0.36e-3),dR_a(4,0.30e-3),dR_a(3,0.22e-3),dR_a(2,0.15e-3),dR_a(1,0.08e-3))
  51. R_2a=set(R_a(11.5,0.82e-3),R_a(11,0.79e-3),R_a(10,0.72e-3),R_a(9,0.64e-3),R_a(8,0.58e-3),R_a(7,0.50e-3),R_a(6,0.44e-3),R_a(5,0.36e-3),R_a(4,0.30e-3),R_a(3,0.22e-3),R_a(2,0.15e-3),R_a(1,0.08e-3))
  52. average(R_2a)
  53. eval(error(R_2a))
  54. count(dR_2a)-count(R_2a)
  55. dR_2b=set(dR_b(11,48.5e-6),dR_b(10,44.0e-6),dR_b(9,40.0e-6),dR_b(8,35.5e-6),dR_b(7,31.0e-6),dR_b(6,26.5e-6),dR_b(5,22.5e-6),dR_b(4,18.0e-6),dR_b(3,13.5e-6),dR_b(2,9.0e-6),dR_b(1,5.0e-6))
  56. R_2b=set(R_b(11,48.5e-6),R_b(10,44.0e-6),R_b(9,40.0e-6),R_b(8,35.5e-6),R_b(7,31.0e-6),R_b(6,26.5e-6),R_b(5,22.5e-6),R_b(4,18.0e-6),R_b(3,13.5e-6),R_b(2,9.0e-6),R_b(1,5.0e-6))
  57. average(R_2b)
  58. eval(error(R_2b))
  59.  
  60.